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Understanding light absorption and scattering mecha¬ 
nisms is of prime importance in optical physics, photon¬ 
ics and plasmonics. Nowadays, the nanofabrication tech¬ 
nologies provide tools to design artifical materials for en¬ 
gineering light-matter interactions. Plasmonic nanopar¬ 
ticles agregates are, for instance, artificial structures 
which have been subject of intense researchs [lj [2]. The 
collective interactions of localized plasmon modes can 
give rise to pronounced Fano resonances on or lead 
to electromagnetic induced transparency [5] or induced 
absorption [6] phenomena. 

In this Rapid Communication, we analyze how strong 
light absorption and scattering can be in arbitrary sets of 
resonant particles when they are cooperatively coupled. 
This problem has been partly considered in the specific 
case of single objects nnn, objects with simple sym¬ 
metry shapes pn ns], periodic arrays of dipoles m or 
dilute distibution of scatterers m ■ Here, we provide a 
general answer to this problem by deriving, in the frame¬ 
work of the linear classical electrodynamics, the funda¬ 
mental upper bounds for light absorption and scatter¬ 
ing for arbitary distributions of resonant particles when 
potentially strong cooperative coupling mechanisms take 
place and when the field radiated by the multipolar or¬ 
ders of particles participate to the coupling. We also 
propose a general inverse strategy to reach these limits. 

Let us start, for pedagogical reasons, by considering 
dipolar systems made with N dipoles of dipolar moment 
Pi spatially distributed at the position in a transparent 
medium of permittivity ey and higlighted by an incident 
field E i nc . The local electric field E i oc measured at any 
point results from the superposition of external incident 
and scattered fields. Therefore it takes the self-consistent 
form duns] 

E ioc (r) = E irac (r) + J 2 /lo^G^r, r,-)p j, (1) 
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is the free space 


Green tensor in the host material defined with the 


unit vector p = p/p , p = r — r', k is the wavec- 


tor in the host material while 1 denotes the unit 
dyadic tensor. By introducing the vectorial fields 
P =t (p( r i), •••,p(r N )), E ioc =* (E toc (ri),...,E Zoc (r N )), 
E mc =* (E inc (ri),...,E inc (r N )) and using the Poynting 
theorem [20], the total power dissipated in this system 
highligted by a monochromatic incident field writes 

EDH22I 


Vabs{u) = ~Im[p\uj).'Ei oc (u))]. (2) 

Using relation 0 and the reciprocity principle [23], this 
expression can be recasted in term of the incident field 

Vabs{u) = -p t .(^y^[/m0]).p + Re[p ] ,i^E inc }. (3) 

where we have set the N x N block matrix of component 
Qij = G 0 (ri, i*j ). Hence, the maximal power than a set of 
dipoles can dissipate under an external highligting reads 
(see Supplemental Material [24] ) 

= ^ E UM _1 -Em C - (4) 

This expression is the first result of the Letter. It repre¬ 
sents the fundamental limit for light absorption by any 
system of dipoles which collectively interact.This expres¬ 
sion depends only on the geometric configuration and on 
the characteristics of incident field. It immediately fol¬ 
lows from it that the maximal absorption cross-section 
for a set of dipoles reads 


max ^ -E jnc /r\ 

abs “4 w v^JE inc | 2 ' W 

This result could seem at first view counter intuitive. In¬ 
deed, with resonant particles we could expect a blowup 
mechanism for light absorption allowing to go beyond this 
limit. However, this intuition is wrong. To be convinced 
of this let us simply examine the case of an isolated par¬ 
ticle of polarizability a = a + ia . In this simple case, 
the dissipated power writes [16] 

V abs C,a?) = ^ | E inc I 2 {a” - Mg I a | 2 }. (6) 
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Then, a straighforward calculation shows that the opti¬ 
mal polarizability is a opt = 3m ^ and it leads to the 
maximal dissipated power P ™ b a *' sln9 {u) = | u2 ™ nh \ 
E inc | 2 which exactly corresponds (using ImQ = ^ 1 ) 
to the value predicted by relation 

Similarly, we can maximize the power 

Vsct = p f .(^y^[/m0]).p (7) 

scattered by a set of dipoles under the constraint V a bs > 0 
to find (see Supplemental Material EH) 

VZT = | 2 _e l c .[imgr\E inc . ( 8 ) 

Comparing this expression to relation Q we see that the 
maximal power than a set of dipoles can scatter is four 
times equal to the maximal power it can absorb. 



Figure 1: Maximal absorption cross-section of a pair of elec¬ 
tric dipoles enlightened (a) in a transversal direction (TE and 
TM polarization) and (b) in a longitudinal direction. In this 
case, the power dissipated by the first dipole is plotted in 
blue dot dashed line while the power dissipated by the second 
dipole is plotted in red dashed line. At subwavelength sepa¬ 
ration distance the first dipole is amplifying. All curves are 
normalized by the maximal cross-section of (resp. the power 
dissipated by) a single dipole, (c) Poynting vector streamlines 
around a dimer of SiC nanoparticles optimized to reach the 
maximal absorption marked by the green disk in Fig. 1(a) for 
a distance d = 7.3 gm and a wavelength A = 10.8 gm. The op¬ 
timal radius is 0.805 gm. The electric field inside the particles 
is two order larger than the unitary incident field. 

The above analyzis gives the ultimate limits that any 
dipolar system can reach given a spatial configuration. 
Surprisingly, these limits are independent of the material 
properties. However, hereafter we show that, in practical 
point of view, these values can be reached by using ap¬ 
propriate nanoparticles. To this end, let us consider the 
general relation p = a between the generalized po¬ 
larizability of dipoles and the regularized version of local 


field defined as E[®®(r») = E inc (r i )+w 2 / u 0 Z) (G o(ri,rj)Pi 

and let us assume, for clarity reasons, that all polariz¬ 
ability tensors are diagonal (i.e. a = diag(a i,...,cqsr) 
with a\ = diag(ai, Xl a^y^a^z)). Then, by inverting this 
relation using the optimal dipolar moments we find for 
i = 1 ,..., N and (3 = x,y,z 

a i8 = _ Popt,i,P __ ( 9 ) 

E inc,i,/3 ~t~ [A o [Co(Uj, Ij ) -P opt,j]p 

j/* 


Here, we must emphasize that, in principle, those opti¬ 
mal polarizabilities do not necessary correspond to polar- 
izabiliies of lossy media so that, generally speaking, the 
optimal absorption can be achieved by combining dipoles 
with gain to lossy dipoles. 

Let us now examine the usefull case of a pair of 
dipoles (separation distance d) along the x-axis: ImQ = 
f Imi G 0 (ri,r 2 ) 

\ /mG 0 (ri,r 2 ) 
a 0 0 
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For an incident field of magnitude Ei nc orthogonal 
to the axis linking the two dipoles, we have for TE 

waves rz x ’ ±TE = fera while for ™ waves 

C’ 1 ™ = This result is P lotted in 

Fig. 1-a. When the separation distance becomes suffi¬ 
ciently large then b ~ a ~ 0 and we see that the power 
absorbed by the pair is twice the power absorbed by 
an isolated dipole. On the contrary, close to the con- 
tact, a ~ g ^[kd - and 6 ~ kd - ^>1] so 

that V2T’ ±TE '™ = V]2 X ' Sin9 - Between these two ex- 
treme regimes the maximum is an oscillating function 
with respect to the separation distance as illustrated in 
Fig.l. in the case of a silicone carbide (SiC) nanoparti¬ 
cles dimer. In this case, the radius of SiC particles has 
been optimized from the optimal polarizability © us¬ 
ing the method described in m to reach the ultimate 
absorption shown by the green disk in Fig. 1(a) at a sepa¬ 
ration distance d = 7.3/mi for a wavelength A = 10.8/ira 
(esic = —2.53 + 0.18z [25]). It clearly appears Fig. 1-c 
that the dissipation of the energy of incident field is due 
to a stong enhancement of field inside the particles. 

When the dimer is lighted along its axis (see Fig. 1- 
b) both dipoles receive a phase shifted incident field. 
















































3 




Figure 2: Maximal absorption cross-section for chains of 
N=10 dipoles enlighted (a) in the transversal direction (TE 
and TM polarization) and (b) in the longitudinal direction, 
(c) Power dissipated per dipole inside optimal regular chain 
with different number N of dipoles separated by a distance 
j = 0.25 for a longitudinal lighting (b). (d) Spatial distribu¬ 
tion of losses and gains inside a chain of 10 dipoles (y = 0.25) 
enlighted along its axis (green point in Fig. 2(b)). The cross- 
sections are normalized by the maximal cross-section of a sin¬ 
gle dipole. The power dissipated in each particle is normalized 
by the total power dissipated inside the chain. 


Let ip = kd be that phase shift and let us denote 
by E inc the magnitude of incident field on the first 
dipole. Then , it is straighforward to see that V abs ’ 7/ = 

I• In tllis C£lSe ’ at cl ° Se se P aration 
distance V^ X,// = \V™ b *' sin9 m This value is the max¬ 
imal power a dimer can dissipate in this configuration. 
At first sight, this result seems to be paradoxal because 
one dipole is in the shadow of the second. But when 
we examinate the optimal losses per particle in near-field 
regime (i.e. ^ < 1) we see Fig. 1-b that one particle dis¬ 
sipates the incident energy while the second not. On the 
contrary, the latter is purely amplifiying. In this regime 
of strong coupling, the higher are the overall losses, the 
higher is the gain in the first particle. This exaltation 
of losses is driven by a cooperative coupling mechanism 
between a lossy and a particle with gain. Even when 
the first particle is lossless (at | ~ 0.195 in Fig. 1-b), 
the power which is dissipated in the second particle and 
therefore the overall dissipated power is more than three 
time larger than the maximal power than a single isolated 
dipole can dissipate. As shown in Figs. 2 similar results 
can be observed with linear chains of dipoles. The spa¬ 
tial distribution of losses inside such chains illuminated 
transversally and along the axis are plotted in Figs. 2-a 
and 2-b in the case of N=10 particles regularly spaced. 
In the longitudinal lighting case, we observe Fig. 2-c an 
almost symmetrical structure that emerges from the op¬ 


timization process with half of the chain which is made 
of dissipative particles and the rest of the chain which is 
made of amplifying particles. As for the variation of the 
maximum absorption cross-section cr^ x with respect to 
the number N of dipoles inside the chain, the numerical 
simulations predict that it increases (Fig. 2-c) linearly 
as a abs X ~ 2 A 3 +1 showing so that apparently there is 
no upper bound for long chains provide had hoc optical 
properties are used. 

In two dimensional networks (Fig. 3-a-and 3-c) we get 
qualitatively the same behaviors than in chains. In con¬ 
trast, in three dimensional networks the optimal power is 
almost independent on the separation distance between 
the dipoles even in near-field regime where the electro¬ 
magnetic couplings are strong. One can speculate that 
both the exaltation and inhibitation mechanims annihi¬ 
late each other. However, the detailed investigation of 
three-dimensional dipolar networks goes far beyond the 
objective of the present work and it is in itsef a problem 
which would necessitate a specific study. 

It clearly appears in Fig. 3, that the absorption cross- 
section of the whole system can, in near-field regime, be 
much larger than the apparent area of domain (shown 
by the green dashed lines) on which the dipoles are dis¬ 
persed. The numerical simulations show that, at sub¬ 
wavelength separation distances, the cross-section can be 
more than one order of magnitude larger than the max¬ 
imal cross-section of a simple dipole.This result outlines 
the importance of collective couplings in this regime and 
show its strong potential for a variety of problem such as 
for instance, the design of subwavelength superabsorbers. 
Nevertheless, it is worth to emphasize that to get an ab¬ 
sorption cross-section which is much larger than the ap¬ 
parent support of dipoles, active media must generally be 
inserted inside the network. When we impose all dipoles 
to be dissipative we observe that the cross-section dras¬ 
tically decreases with the separation distance to become 
comparable with the maximal cross-section of a single 
dipole (shown by the horizontal blue dashed lined in Fig.3 
) near the contact. 

So far, we have only considered couplings between par¬ 
ticles described by simple dipoles. Hereafter we consider 
the most general situation where higher order modes (ie. 
multipoles) are taken into account [26]. The electromag- 
ntic field inside a medium of refractive index rih can be 
expressed in term of ingoing (-) and outgoing (+) vector 
spherical wave functions (which form a complete basis) 

*-(£) <l0, 

where we have adopted the usual convention for the mul¬ 
tipolar index (m, n) which are replaced by a single index 
p = n(n + 1) + m and where q set the polarization state 
(i.e. q = 1 for TE waves and q — 2 for TM waves). The 
outgoing wave functions are solutions of Maxwell’s 














4 




Figure 3: Maximal absorption cross-section of two- 

dimensional (resp. three dimensional) random networks of 
20 dipoles and (resp. 40 dipoles) for a perpendicular and a 
longitudinal highligting with respect to the surfacic fraction in 
dipole.Here the results of 50 realizations (with a uniform prob¬ 
ability density) are plotted for purely absorbing dipoles (red 
curves) and arbitrary dipoles (black curves).The blue dashed 
curve shows the maximal absorption cross-section for a single 
dipole. The green dashed curves denote the apparent area of 
domain in which the dipoles are distributed. 


equation (using the e lujt convention) 


V x E+ - iujfiH.tr 


U s 

pq 

V x H+ = -iue E+ + E® 


pq 


pq 


( 11 ) 


with the source term ^ = 


rpS' 

E s? 

pq 


H 


. Now, let us intro¬ 


duce the following real spherical harmonic functions 


X 


± 

pq 


' 7f(^W + ^pq) if TO<0 
< 'ippq if m = 0 

.72 (^-^) if m>0 


( 12 ) 


and let us define the regular harmonics ( pq = V '"T'''' . 
By definition, the incident field can be decomposed, in 
the complete basis of real spherical harmonic functions, 
as 


V ; mc( r ) — y^ a fnc.iCpg( r r i)- (13) 

pq 

As for the local field it takes the form 


^Zoc(r) = Ipinci r) + 5Z[ a i 9 Cpq ( r - r i) + A TXpq (r - r i)j- 
pq 

(14) 

where the second term of rhs is the total scattered field 
which corresponds to incoming field scattered by all par¬ 
ticles in direction of the i th scatterer and the outgoing 


field diffracted out this particle. Since the incoming field 
is related to all outcoming fields by a linear relation of 

the type a^ q = ^ Tf / ,P q A^ q , where the T pq,p q are 

j^ip'q' 

the components of a translation operator (a propagator) 
m between the i th and the j th particle, the local field 
takes the form 

V’Zoc(r) = 

£X+ (r - rOdaSL* + | + ^f) 

pq j^ip'q' 

+Ex p -,(r - *)(§«&,« + \ £ ET?/’ pV ). 

pq j^ip'q' 

(15) 

Then, using the orthogonality relations for the functions 
Xpg, the power dissipated in each particle can be rewrit¬ 
ten as the net flux 

Pabs,i = El" I hind + \ E E A P '«' + Af | 2 

pq j^ip'q' 

+ lKn\ i + lEE^r v Af' I 2 } 

j^tp'q' 

(16) 

accross a surface surrounding the particle. By develop- 
ping this expression we get 

Va bs , = -E(l A i 9 I 2 +Re[AT\aZ c ,) 

+EE^ 9 ,pV T V ]}- (1?) 

j^ip'q' 

Summing the losses over all dissipating objects we get, 
after a straightforward calculation, the total power dissi¬ 
pated by the system 

'Pabsi^) = Y^'Pabs,i = —AFV .A + Re[ A^.B] (]_g) 

and the power scattered by it 

Vscat(u) = At.V.A. (19) 

In these two expressions, A and B are, in a system of N 
scatterers, the block vectors defined with the subvectors 
A i and = -(a^ c i ) while V = 1 + u + ut is the block 
matrix defined with the sub-blocks elements U^ = T^. 
Following the same reasoning as used for a set of dipoles 
we obtain the upper bounds for the dissipated and the 
scattered powers inside any arbitrary system of particles 
which are collectively coupled (see Supplemental Mate¬ 
rial [24]). 

These results are illustrated in Fig. 4 in the specific 
case of a dimer of nanoparticles. We note that, similarly 
to the results obtained for single objects [8], the maximal 
cross-section is an increasing function with the number of 
multipolar orders. In fact the number of resonant modes 
of the whole system increases with the multipoles order 
creating so new channels for dissipating or scatter light. 

In conclusion, we have analytically derived the ulti¬ 
mate limits for light absorption and sctattering by a sys¬ 
tem of N point resonant multipoles for a given geometry 
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Figure 4: Maximal absorption cross-section for a dimer 
of nanoparticles versus the separation distance.The green 
dashed line corresponds to the dipolar (electric) case while 
the blue, purple and black lines stand for the contribution of 
multipolar orders one (dipolar electric+magnetic), two and 
three, respectively. 


optimal problem. Find X opt , 


'PabsO^opt) Max{P abs (X)} 

With Pabsfi^opt') ^ 0 


( 20 ) 


where, according to relation (3) the absorbed power 
writes under the form V a bs — — A"^.V.A + i?e[A^.B] with 
A = p, V = ^-^-[ImQ] and B = inc in the dipolar 
case. By decomposing this expression with respect to 
the real and imaginary parts of those vectors we have 
equivalently 



of their spatial distribution. We have demonstrated that 
these limits are independent of the optical properties of 
materials but they can be reached only by combining 
lossy with active media. Beside the mechanisms of cou¬ 
pling between light and single objects, the cooperative 
interactions mechanisms in systems of coupled dissipat¬ 
ing and amplifying resonant particles offer a supplemen¬ 
tary degree of freedom to tailor light-matter interactions. 
Those mechanisms pave the way for promising scientific 
and practical applications in optics. 

Supplemental material 

In this supplemental document we give supplemen¬ 
tary informations concerning the derivation of the op¬ 
timal power absorbed or scattered by a set of optical 
resonators. The derivation is general and can be used 
either in the dipolaror the multipolar case. 

I. Maximal absorption 

The maximal power which can be absorbed by such a 
system can be found by solving the following constrained 


Figure 5: Illustration of two optimization problems.On the 
bottom, the paraboloid corresponds to all physically admissi¬ 
ble dissipated powers inside a system over the space of con¬ 
figuration (brown plane). The orange area on the upper part 
of this praboloid surface corresponds to lossy systems while 
the lower part correspond to amplifyiing systems. The blue 
sphere on the top of this surface denotes the maximal dissi¬ 
pated power while its orthogonal projection on the configura¬ 
tion plane stands by a blue cross marks the optimal general¬ 
ized polarizability. The upper part of this mapping illustrates 
the maximization of the scattered power. The green paraboid 
stands for all physically admissible scattered powers. The 
optimal generalized polarizability under the constraint of a 
positive dissipated power is shown by the red sphere. 


By construction A is a symmetric matrix so that ex¬ 
pression (21) is a quadratic polynomial. Hence, the 
maximal dissipated power can be obtained by solv¬ 
ing constrained optimal problem (1) using the classical 
Kuhn-Tucker (KT) conditions (see S. Boyd and L. Van- 
denberghe, Convex Optimization (Cambridge University 
Press, Cambridge, 2004). In the present case, this is 
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equivalent to cancel the gradient <r ^’‘ = —A.X + B so 
that 


X opt = A~ 1 .B 


1/V- 1 0 

2 \ 0 V- 1 


i?e[B] \ 
Im [B] ) 


( 22 ) 


Here below we give another demonstration for the op¬ 
timal powers absorbed or scattered by a N-body system 
following a purely algebraic approach. By construction V 
is a postive definite real matrix. Thus the two following 
scalars 


It follows from the definition of vector X that the optimal 
generalized dipole vector field reads 


Qi = W - Vv-^.V-tA. - Iv^B) (30) 


A opt = Iv^B. (23) 

By inserting this expression into relation (3) of the Letter 
we get the maximal absorption 

= Ifit.V-^B (24) 

which is equivalent, accordng to the definition of different 
vectors and tensor fields to 

PSTM = Einc , (25 ) 

II.Maximal scattering 

The maximal power which can be scattered by an ar¬ 
bitrary system can be found by solving the following con¬ 
strained optimal problem. Find X opt , 

f Vset (3^opt ) ^ Max{V sct (X)} 

< x (26) 

[with Vabsi^-opt) > 0 

where V ac t = AhV.A. Using the matrix A and the 
vector B introduced above the scatterd power writes 
equivalently V sc t = ^X^.AX . Then, to solve the 
constrained optimal problem (6) we introduce the La- 
grangian £(X, A) = |XlAX + W a bs{ X) where A > 0 
denotes a positive KT-multiplier [? ]. Thus, according 
to the KT conditions, the optimal solution must can¬ 
cel the gradient £x = AX + A(£> — AX) of this La- 
grangian. This condition implies that X = jzr[A~ 1 B. 
Beside, the constraint V a bs — 0 gives the optimal KT 
multiplier A = 2 so that case X opt = 2A~ l B. It follows 
from the definition of vector X that the optimal dipole 
vector field reads 


Aopt = V -1 .B. (27) 

Finally, by inserting this expression into the relation (6) 
of the Letter we get 

v ma X _ V -1 .B (28) 

which is equivalent to 

VT c 1 x = ^_E \ nc .[Img]-\V inc . (29) 

III.A simple alternative for the derivation of 

optimum 


and 


0 2 = (A t .-B t V- 1 ).V.(A.-V- 1 B) (31) 

are positive. From these expressions the absorbed and 
the scattered power given in Eqs.(3) and (6) of the Letter 
write 


Tabs = -Bt.V-i.B _ Ql 


(32) 


and 


Vscat — B^.V 1 .B — Q 2 — 2 Vabs: 


(33) 


respectively. According to (32) V a bs is bounded by 


'T)max 
’ abs 


= iBt.V 1 


B and from the definition of Qi we see 
that this value is reach when A 0 jrt = |V 1 .B. Following 
a similar reasoning for the scattered power we have from 
@ VIZ'f < BFV FB. Then, by choosing A = V : B 


we see, according to expression (31), that Q 2 = 0 and 


Vabs = 0 and we reach the upper limit for the scattered 
power. 
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